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g ■ Abstract 
> 

^ ■ We re-examine the very stringent hmits on the axion mass based on the 



strength and duration of the neutrino signal from SN 1987A, in the light of 
new measurements of the axial- vector coupling strength of nucleons, possi- 
ble suppression of axion emission due to many-body effects, and additional 
emission processes involving pions. The suppression of axion emission due 
to nucleon spin fluctuations induced by many-body effects degrades previous 
limits by a factor of about 2. Emission processes involving thermal pions 
can strengthen the limits by a factor of 3-4 within a perturbative treatment 
that neglects saturation of nucleon spin fluctuations. Inclusion of saturation 
effects, however, tends to make the limits less dependent on pion abundances. 
The resulting axion mass limit also depends on the precise couplings of the 
axion and ranges from 0.5 x 10~^ eV to 6 x 10~^ eV. 
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I. INTRODUCTION 



Peccei-Quinn (PQ) symmetry ||T]| continues to be an attractive solution to the strong CP 
problem. However, theory and laboratory experiment give little guidance on the PQ symme- 
try-breaking scale, /a, and therefore on the mass of the pseudo-Goldstone boson associated 
with PQ symmetry, the axion [^: 

rua ~ 0.62 eV ■ 10^ GeV//a . (1) 

Astrophysical and cosmological arguments have been very powerful in excluding values of 
the axion mass, allowing only the mass range from about 10~^eV to about 10~^eV f^. 

Axions in this mass range would have been produced cosmologically by a coherent, non- 
thermal mechanism Because of this, axions in this mass range would be cosmologically 
significant, have small velocities, and behave as cold dark matter (CDM), in spite of their 
small mass. Depending upon cosmological and particle physics parameters, axions could 
contribute to the closure mass density today for nia — 10~^eV to 10~^eV [|]. We recall 
also that the CDM scenario for structure formation is at present the most promising 
with axions and neutralinos being the leading CDM candidates. If axions provide the bulk 
of the dark matter, they must comprise a significant fraction of the dark halo of our own 
galaxy and a large-scale experiment is underway to detect halo axions of mass 10~^eV 
to 10-5 eV [|]. 

The most stringent astrophysical bound on the axion mass is that derived from obser- 
vations of neutrinos from SN 1987A, excluding masses greater than about 10~^eV P,p!o|. 
It is based upon the detection of 19 neutrinos from SN 1987A by the Kamiokande II (KII) 



TT| and Irvine-Michigan-Brookhaven (1MB) [|T^] water- Cherenkov detectors. According to 
the standard model of Type II supernovae, namely the core collapse of a massive star, these 
neutrinos were emitted during the early cooling phase of the nascent neutron star associated 
with the appearance of SN 1987A. Indeed, the observed neutrino flux and energy spectrum 



are consistent with this picture [13|. The emission of axions would have hastened the cool- 
ing process, leading to fewer events over a shorter time. The axion's couplings to ordinary 
matter are proportional to the axion mass, so consistency with the detected neutrino burst 
leads to an upper limit to the axion mass, which is estimated to be around 10~^ eV. 

A number of questions have been raised about the mass limit based upon SN 1987A. They 
include the possible suppression of axion emission by many-body effects which are likely to 
be important in the deep interior of a neutron star [p!^ -[l6[], enhancement of axion emission 



due to the possible presence of kaons and thermal pions at the centers of neutron stars |]T7 



and the implications of recent measurements of the strengths of axial-current couplings to 
nucleons, which may lead to a significant suppression of the axion-neutron coupling relative 
to estimates based on the naive quark model (NQM) |T6[] . 

The purpose of this paper is to address these issues; its outline is as follows. Sect. II 
is devoted to a discussion of the axion-nucleon couplings and a description of the input 
microphysics of the numerical protoneutron star models. The axion emission rates are 
reviewed in Sect. Ill, including many-body effects, with particular attention to the possible 
damping effects of nucleon spin fluctuations. Sect. IV discusses cooling models and the 
impact of axion emission on the theoretical predictions for the neutrino bursts detected by 
KII and 1MB. We finish with a summary of our results in Sect. V. 
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II. INPUT PHYSICS 



A. Axion Couplings 

The most important axion-emission process in a hot, young neutron star is axion brems- 
strahlung when two nucleons colhde. The rate for this process depends upon the axion's 
couphng to nucleons. The relevant part of the interaction Lagrangian is 

>Cint = TT^V'iTMTsV'i , (2) 

where the index i denotes a neutron or proton, and we can write the axion-nucleon couplings 
Qai in the form 

Qai = CiUli/^fpQ/N) = CiUli/fa . (3) 

is the color anomaly of the PQ symmetry. The dimensionless couplings Ci of Eq. are of 
order unity, and related by generalized Goldberger-Treiman relations to nucleon axial- vector 
current matrix elements via the PQ charges Xg of the light quarks q = u,d, s: 

Cp = [xu - 1/(1 + Z + w)]Au + [xd - z/{l + Z + w)]Ad + [x, - w/{l + z + w)]As; 

Cn = [xu - 1/(1 + z + w)]Ad+ [xd - z/{l + z + w)]Au + [xs - w/{l + z + w)]As; (4) 

where we use for the light-quark mass ratios 

z = rriu/md — 0.565 , w = mu/ms ~ 0.029 . (5) 



The PQ charges of the quarks are model dependent: they vanish for the KVSZ axion [18 



Xu = Xd = Xs = 0, whilst for the DFSZ axion [|19| they can be written in terms of an angle 
f3 which is related to a ratio of Higgs vacuum expectation values: 

Xu = sin^/3/3, Xd = Xg = cos^/5/3 (6) 

The Aq in Eq. (^ quantify the axial-vector current couplings to the proton: 

Ag5^ = (p|g7^75g|p) (7) 

where S*^ is the proton spin. Similar expressions hold for the neutron, with matrix elements 
related by an isospin reflection: Am„ = Ad, Adn = Au. 

The Aq are non-perturbative quantities whose values must be taken from experiment. 
Neutron (3 decay and isospin invariance constrain An — Ad = (7^ ^ 1.25, whilst hyperon (3 
decays and flavor SU{2)) symmetry for the baryon octet yield Am + Ad — 2As ~ 0.682. The 
best determinations of the third combination of the Ag are obtained from spin-dependent 
deep-inelastic electron and muon scattering off nucleons. Recent analyses give [EO]: 



Au= 0.80 ±0.04 ±0.04 
Ad= -0.46 ±0.04 ±0.04 
As = -0.12 ±0.04 ±0.04 



(8) 



3 



where the first error is a statistical error, and the second is an estimated systematic error. 
The values in Eq. (§) are quite different from those estimated in the NQM, which does not 
estimate axial-current matrix elements reliably. As can be seen in Eq. the measurements 
of Eq. (^) and the PQ charges of the light quarks determine the axion-nucleon couplings. 
Values for the KSVZ and DFSZ axions are shown in Fig. |l|. 

It is perhaps worth warning the reader at this point that the above values of the Ag have 
been obtained for nucleons that are free, or in light nuclei. It is possible that the axial- vector 
current couplings may be different for nucleons in a dense medium. However, large effects 
on these couplings are not seen in conventional nuclear physics, and we believe that other 



many-body effects on the axion emission rate are dominant, as we discuss in Sect. Ill C 



B. Protoneutron Star Physics 



In order to obtain limits on the axion production in SN 1987A and thus on the axion 
mass, we performed numerical simulations of the evolution of newly formed neutron stars. 
The simulations started a few milliseconds after core bounce and supernova shock formation 
and followed the Kelvin-Helmholtz cooling by neutrino emission until the lepton-rich and hot 
protoneutron star had evolved to the final cold and neutronized state several ten seconds 
later. Cooling sequences were computed with and without axion emission included, and 
the corresponding neutrino luminosities and spectra were used to derive predictions for the 
associated KII and 1MB detector signals. 

The initial models of the protoneutron star were constructed with profiles of the electron 
concentration Ye = Ue/ub {fie is the electron number density, Ub the baryon number density) 
and of the temperature T that were very close to those obtained in detailed hydro dynamical 
calculations of stellar core collapse [^. The collapsed stellar core right after core bounce 



is gravitationally only weakly bound. This is expressed by the fact that the ratio of the 
gravitational mass at the beginning of the simulation, M^^^g, to the baryonic mass, Mns b, is 
only slightly less than unity. We scaled the temperature profiles from given numerical models 
of post-collapse cores such that M]^^^^/ M^s,h = 0.97 in all of our initial setups. The baryonic 
mass of the neutron star depends on the mass of the progenitor star of the supernova. 
Simulations with varied protoneutron star mass [0,^] showed that the neutrino signal 
from SN 1987A can be best reproduced by models with a mass of Mns,b — 1-5... 1.6 M©. 
The reference model of the simulations presented here was therefore chosen to have Mns,b = 
1.53 Mq and an initial gravitational mass of M^^ g = 1.49 Mq. The electron number fraction 
had a maximum at the center of the star where Y^^ ~ 0.29. The central temperature and 
density at the beginning of the simulations were T^. = 24MeV and p]. = m x 0.485 fm~^ = 
8 X lO^'^gcm"^ ~ 3.25pnuc5 respectively, with pnuc being the nuclear matter density and m 
the common nucleon mass. 



The computations were done with a general relativistic stellar evolution code [22 



em- 



ploying equations of state developed for the case T = by Glendenning ||23| and extended 
to finite temperatures by including thermal corrections to the energy density and pressure 
of both stellar plasma and neutrinos Two different equations of state were employed. 
As the standard case, we used a hyperon equation of state (case 2 in EOS B in p2[ ) 
which takes into account n, p, e^, and fi~ in the nuclear medium, and hyperons (A, S...) 
and A-resonances as additional hadronic states at densities above about 2pnuc- Due to the 
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formation of these additional baryonic degrees of freedom, the hyperon equation of state is 
"softer" at very high densities than a non-hyperonic equation of state involving only n, p, 
and fi~ . An equation of state of the latter category (case 5 in [Q, EOS A in |Q) was 
used for comparative computations in the work presented here. 

The equations of state employed here do not predict the occurrence of pions or pion con- 
densates in the supernova core. The possible importance of pionic excitations was pointed 
out in a paper by Mayle, Tavani, and Wilson [^. The presence of a significant number of 



negative thermal pions influences the axion emission from forming neutron stars by produc- 
ing axions via pion-axion conversion processes tt~ + p ^ n + a (Sect. [1IIB| and jl^). In order 



to estimate the impact of the latter on the protoneutron star cooling, we computed a number 
of models using the rather crude assumption that one tt~ per nucleon is present at densities 
beyond about two times nuclear matter density. The local abundance Y^^- = n^-/nb was 
simply prescribed by F^- = Yn ■ {1 + tanh [2(p — Pnuc)/Pnuc]} /2 with pnuc = m x 0.15 fm~^. 
This prescription is ad hoc and effects of pions on the nuclear equation of state were not 
taken into account at all. Also, we ignored the temperature dependence of the abundance 



of thermal pions. Our choice of Y^- of order unity was motivated by the results of Ref. pi 
where an abundance of Y^^- < 0.6 (and a total pion abundance of Y^^ ^ 0.8) was obtained for 
the conditions of temperature and density in a neutron star model at 3 s after core bounce. 
The possibility of a pion condensate, however, was excluded for the supernova conditions and 
the pion dispersion relation assumed in Ref. Since our whole treatment of the pion case 
was very approximate only, we also did not make any attempt to include the contributions 
of pions to the neutrino scattering opacity of the nuclear medium (see also [|4 



As for the neutrino transport, equilibrium diffusion was assumed for all types of neutrinos 
(z/e, z/e, z/^, z/^, Ur, ^r) which is & good approximation, because the hot matter of the protoneu- 
tron star is very opaque against neutrinos by neutral-current neutrino-nucleon scatterings, 
1/ + N ^ u + N, and charged-current /3-reactions, Ue + n ^ e~ + p and z/g + p — >■ e"*" + n. For 
details of the technical implementation and the "standard" description of neutrino-matter 
interactions, see . A possible reduction of the neutrino opacity by many-body correlation 
effects, e.g., a suppression of the uN scattering cross section by rapid spin fluctuations due 
to frequent nucleon-nucleon collisions ( [|T^ - [l6| , ^ , p6| and Sect. |III C|) , was included only in a 
comparative model run to reveal the change of the axion mass limit. We want to emphasize 
here that the central density of our reference neutron star model with M^s,h = 1.53 Mq is 
always less than 8.8 x lO^^gcm"^. Only in a relatively small, central part with a mass 
< 0.5 Mq does the density become higher than about 2pnuc and hyperonization sets in. But 
even at the center of the star the hyperon abundance Ihyp = '"-hyp/'^.b never rises above 
Yhyp — 0.25 (see Fig. 7 in |2^). Therefore we consider the disregard of modifications of 
the neutrino opacity due to the presence of hyperons as acceptable for the models dis- 
cussed in this work, and the assumption that the neutrino opacity is produced mainly by 
interactions of neutrinos with n and p should yield a sufficiently accurate description. 

In order to evaluate our protoneutron star cooling models for the prediced neutrino 
signals in the KII and 1MB detectors, we folded the computed spectral Ue number flux 
with the detector efficiency functions and the cross section for z/g absorption on protons. A 
distance to SN 1987A of D = 50 kpc was assumed. For details of the evaluation procedure. 



see Appendix C of [22 
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III. AXION EMISSION RATES 



A. Nucleon-Nucleon Axion Bremsstrahlung 



As long as the widths of the nucleoli states are small compared to the temperature, one 
can evaluate the matrix elements for nucleon-pair bremsstrahlung of axions by using free nu- 
cleon states. A one-pion exchange (OPE) potential is likely to be an adequate starting-point 
for describing the two-nucleon interaction. As a function of the four-momentum transfer be- 
tween the nucleons, k = {kg, k), it can be written as 



OPE 



;k, 0-1,0-2) = - — 




(0-1 ■ k) (0-2 ■ k) 



fc2 



[Ti ■ T2J 



(9) 



Here, / ~ 1 is the pion-nucleon coupling constant. 



is the pion mass, and o-,- and r,- 



are spin and isospin operators for the two nucleons, respectively (j = 1,2). The resulting 



30 



matrix element A4 was first calculated in Ref. for degenerate nucleons and later in 
for arbitrary nucleon degeneracy. The lowest-order energy-loss rate per unit volume due to 
axion emission, is then given by the phase-space integral 



(1) 



2c.(2vr)3 2E,(2vr)3 ''^^^^^^ " ^'^^^ " 
S J2 1^ r(2vr) V(pi +P2-P3-P4- ka) , 



(10) 



spins 



where pj = {Ej,pj) are the four-momenta of the initial-state {j = 1,2) and final-state 
{j = 3,4) nucleons, and ka = (ci;, k^) is the axion four-momentum. Furthermore, fj is the 
nucleon occupation number in state pj (j = 1, ■ ■ ■ , 4), and S is the usual symmetry factor. 

In Ref. the 15- dimensional integration in Eq. (^) was performed exploiting the fact 
that varies only slightly in the range where the integrand contributes most. Neglecting 
this variation induces an error of less than a factor 2. Introducing the thermal average 





(P2 - 


-P4) 


(P2 - 


- P3 


( 


|P2 


- P4I 


|P2 - 


PsI 



(11) 



which can be shown to take the values and 1.0845 in the limits of degenerate and non- 
degenerate nucleons, respectively, the result can be written as 



g« = 64 



/ 



^2.5^6.5 



;i - i/^)9lj{yn, Vn) + (1 - i/^)9lpl{yp^ Vp) 



4(15-20 fgL + g. 



ap 
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iiyn,yp) + 



4(6 - 40 fgan + gap 



9 



;i2) 



Here, yi is the dimensionless non-relativistic version of the nucleon chemical potential fif. 
Hi = {f^i — ^i)/T, and the dimensionless function I{yi,y2) can be fitted to within 25% 
accuracy by the analytic expression |^| 
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2.39 X 10^ i^e-y'-y' + 0.256-2^1 + 0.25e-y^) + 1.73 x 10^(1 + \y\) 



-6.92 X 10^(1 + \y\y^^^ + 1.73 x 10^(1 + \y\y 



with y={y^ + y2)/2. 



-5/2 



(13) 



B. Pion-Axion Conversion 



If pions or kaons are present in a supernova core, additional processes such as 7r~ + p ^ 
n + a can contribute to axion emission. Since it is uncertain whether a pion condensate 
can form in the hot protoneutron star we consider only thermal pions here. The 

corresponding lowest-order perturbative energy emission rate per unit volume was found to 
be 



n„-n. 



V 1 



(14) 



where QaN is a combination of axion-proton and axion-neutron couplings, 



-2 
9aN 



\{g 



I 2 
9 an 



■ 9ap9a 



(15) 



and n^- and Up are the number densities of vr and protons. If the pion abundance is 
comparable to the nucleon abundance, energy loss by pion-axion conversion will dominate 



over nucleon-pair bremsstrahlung of axions by more than a factor of 10 in the pertur- 
bative approximation, i.e., when saturation due to many-body effects in the dense medium 
(see Sect. Ill (J| ) are ignored. We have considered the axion mass bounds resulting from 
the lowest-order pion-axion emission rate as well. Due to the clear dominance of the axion 
emission from pion conversion processes in the perturbative approximation, we neglected 
the production of axions by nucleon-nucleon bremsstrahlung when 7i~ were assumed to be 
present in the cooling neutron star. 



C. Many-Body Effects 



Recently the possibility was discussed |jTj-|T^,Q that near nuclear densities many-body 
effects might suppress substantially the actual energy-loss rate in axions compared to the 
lowest-order result Eq. (|T2|). A suppression is likely to occur also for the rate of Eq. (p!^. In 
order to discuss this issue and to compare different approaches, it is convenient to reformulate 
the loss rate in terms of the structure function formalism adopted in these references. In the 
limit of non-relativistic nucleons, axions couple exclusively to the dynamical spin-density 
structure function (SSF), as can be seen from Eq. (|^). Its long- wavelength limit is defined 
as 



4 r+oo 

SJuj) = lim / c/te*^* {ait, k) ■ cr(0, -k)) . 

k^o Srib J-oo 



(16) 
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Here, cr(t,k) = J d'^re~^^'^cr(t,r) (with V the normahzation volume) is the Fourier 
transform of the local nucleon spin-density operator <T(t,r), nj, is the baryon density, and 
{uj, k) is the four-momentum transfer to the medium. The expectation value (■ ■ ■) in Eq. (|16|) 
is taken over a thermal ensemble and the states involved are normalized to V. For a single 
species of nucleons whose coupling to axions is given by Ci [see Eqs. (0)-(§)] Qa can then 
be written as [p!^ , p6[ 

(17) 

For two nucleon species one can absorb the constant Cf into the definition of by 
multiplying cr(t, k) in Eq. (0) by [1 + r^it, k)] Cp/2 + [1 - T^it, k)] CJ2. Here, T^it, k) is 
the third component of the isospin operator T{t, k) which is defined analogously to cr{t, k). 
For the qualitative discussion in the rest of this section, it is sufficient to focus on the case 
of a single nucleon species if not stated otherwise. 

It is obvious from Eq. ([T6|) that only interactions which do not conserve the total nucleon 
spin can lead to non- vanishing values of 5*0- (cj) at uj ^ 0. As can be seen easily, the OPE 
potential Eq. has this property. For the following, it is convenient to introduce the 
lowest-order effective spin fluctuation rate, 

il/2 




m^rib {v) ~ 32 MeVpuTl/,^ , (18) 

where (v) is the average relative velocity between two nucleons, pi4 is the mass density in 
units of lO^^gcm"^, and Tio = T/lOMeV. In the limit of non-degenerate nucleons, the 



lowest-order contribution to the SSF takes the form 15 



Si'\u;) = ^so{uj/T), (19) 



where K = VItx^^/"^ {T /m)^^'^ / (v) ~ 2.7 Q The dimensionless bounded function sq{u/T) 
has been given in Ref . |T5[ . Adopting Eq. ( pj] ) in Eq. (p!7| ) leads to a rate Qa which, for the 
case of two nucleon species, coincides with Eq. (|T^ in the non-degenerate limit, i/n, Up ^ — 1- 
Note that, according to Eq. (0), I{yi,y2) oc e^^^^^ oc nin2T~^ in this limit, where rii is the 
number density of species i. 

However, Eq. (|TUp is unphysical in the limit of both small and large energy transfers, 
which can be seen as follows. First, one can derive the sum rule, 

'^'^5.(^) = i+-^(5:a,-a,) , (20) 



oo 



2ti 




where the sum over all nucleon pairs accounts for possible spin correlations among different 
nucleons. At least at low densities, these correlations can be neglected compared to the first 



^The relation between our r^^'' defined in Eq. ( p^ ) and the quantity Ta used in Ref. [|T^] is 
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term. In any case, the finiteness of the integral in Eq. clearly shows that the infrared 
singularity in Eq. ([T9|) is unphysical. In fact, higher-order effects are expected to regularize 



this singularity at low energy transfers. As a first qualitative guess, it was suggested |15| to 
substitute cj"^ by (cj^ + a^F^)"^, where a is a dimensionless number of order unity which 
can be chosen to satisfy the sum rule Eq. (pOf). 

Secondly, one of us (G.S.) recently derived the analog of what is usually called the 
f sum rule, 

where Ht = | Vij is the "tensor component" of the nucleon interaction Hamiltonian, 
defined as the tensor component of the two-nucleon interaction potential summed 



over all pairs. It is this component which violates local nucleon spin conservation and, 
according to Eq. (^I]) , determines the width of the SSF via the effective spin fluctuation rate 

r^^) defined in Eq. ([T8|) can be shown to be the first-order approximation to To- in a dilute 
medium. It turns out that in the case of two nucleon species the f sum, Eq. (|2T|) , diverges 
when Eq. (^) is substituted for the SSF [^. This can be traced back to the unphysical 
behavior of the dipole-like OPE potential Eq. (P) at small distances. More realistic potentials 
which account for hard-core repulsion lead to a fall off of at high a; which is stronger 
than C(j~^, thus assuring f sum integrability. 

A plausible modification of the SSF is given by 



SAu;) = -jJ^s{u;/T), (23) 



where the continuous and bounded function siu /T) is even and satisfies s(0) = 1. This 
expression has the right limiting behavior in the classical regime which obtains for uj <^ 
T For oj < Fcr, multiple collisions become important and lead to a suppression of 

Scr{uj) compared to the lowest-order approximation Eq. (|19D, which is known as the Landau- 



Pomeranchuk-Migdal effect ||33| , |3^ and ensures integrability of the sum rule Eq. (|20|) . In the 
classical bremsstrahlung limit of hard collisions one would have s{x = uj/T) = 1. Quantum 
corrections require that s{x) —>■ for x —>■ oo sufficiently fast, as can be seen from the f 
sum rule Eq. (^). The exact shape of s{x) depends on the nucleon interaction potential. 
In particular, the high x behavior is governed by the small-distance regime which is usually 
dominated by a hard-core repulsion. 

Modifications of s{uj/T) at a; > T do not have a big influence on the loss rate Eq. (p!7D, 
because of the exponential factor. The finite-width modification is also unimportant as long 
as Fg. < T, i.e., in the dilute medium. In this regime, we have Fq. ^ F^^^ and Qa — Q'^^ oc 
UhT^^T^ [see Eq. (]T^)]. For ^ T, in contrast, Qa/inbT'^) would start to decrease with 
increasing Fo- This would be the case if Fg. ~ F^^-* up to the highest densities [see 



Eq. (|T8|)]. However, in Ref. [^U] it was argued that Fg- is likely to saturate at some maximum 
value F™'^^ ^150 MeV. This also implies that neutrino opacities are suppressed by less than 
~ 50% compared to the lowest-order opacities, ensuring consistency of the observed and 
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simulated cooling time scales for SN 1987A |2§. For < 150 MeV, Qa ^ Qi^^ IS a good 
approximation. For a first improvement to account for saturation of Fo-, we therefore set 



with F;;'^" = 2TrW (24) 

in the numerical simulations, with taken from Eq. ([12|). The maximum value of the spin 
fluctuation rate, F™^"", is coined in terms of the average interaction energy W of a. nucleon 
in the nuclear medium. This definition is motivated by Eq. (p^ and implies the relation 
W = -{HT)/{nhV). 

We have made no systematic attempt to include the effects of a saturation of nucleon 
spin fluctuations due to many-body interactions in case of the pion-axion conversion process. 
However, irrespective of whether the process involves real pions or virtual pions from a 
bystander nucleon, axions are emitted by the fluctuating nucleon spin whose fluctuation 
rate F^- is expected to saturate according to T^^"^ of Eq. (|^). Inclusion of saturation effects 
is therefore likely to leave axion bounds considerably less sensitive to the abundance of pions 
than suggested by the lowest-order energy- loss rate of Eq. (|1^. 



IV. AXIONS AND PROTONEUTRON STAR COOLING 



A. Models 

In this section we discuss our protoneutron star cooling calculations which were per- 
formed with a systematic variation of the axion-nucleon couplings and with different input 
physics in the neutron star modeling. The considered axion-nucleon coupling constants Qai 
are of the order of 10~^° and thus axion opacity plays no role ||35|| . 

In our first sequence of cooling models, the loss of energy by nucleon-pair axion brems- 
strahlung [Eqs. ([ID and (|13D] was investigated in dependence of the two parameters gap 
and X = Qan/ 9ap which were chosen from the intervals < S'ap/lO"^" < 10 and —3 < x < 
3, respectively. Suppression or saturation of the axion emission rates was not taken into 
account. The nucleon degeneracy parameter ^ of Eq. (|lT]) was set to 0.5, but comparative 
calculations with ^ = 1.0845 (nondegenerate nucleons) and ^ = (very degenerate nucleons) 
revealed only a very weak dependence of the results on the particular value of S,- 

In a second sequence of models we repeated these cooling calculations with many-body 
effects taken into account according to Eq. (^) where Q^^^ is from Eq. ([T^p and F^^^ is defined 
in Eq. ([T8|). The average interaction energy per nucleon, W , was chosen toheW = 10 MeV, 
corresponding to F™**^ ~ 60 MeV. A finite value of W leads to a saturation of the axion 
emission rate for large values of F^^\ i.e., at high densities and/or high temperatures. Our 
"standard" case without saturation is formally recovered for F™*^^ <x. W ^ oo. We also 
performed comparative computations for the values W = b MeV and ly = 20 MeV. 

In a third set of models we considered the effects of the presence of a large number of 
negative pions in the nuclear medium. Pions were included in an ad hoc manner as described 
in Sect. |11 lj| by assuming 0.5 < Y^^^ < 1 for the number of tc~ per nucleon at densities above 
Pnuc- The rate of energy loss via the pion-axion conversion reaction is given by Eq. (|n|). 



Qa = gi') Min 
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Because of the 10-50 times larger rate compared to nucleon-nucleon axion bremsstrahlung, 
a range of coupling parameters < ^a7v/10~^° < 2 was explored. 

In addition to these three sets of models and the variations of ^ and W mentioned 
above, we replaced the hyperon equation of state (EOS B) by an ordinary n-p equation 
of state (EOS A) in order to reveal the differences for the axion production during the 
Kelvin- Helmholtz cooling of the protoneutron star when hyperons and A-resonances are not 
present at densities beyond 2pnuc- Furthermore, the influence of a reduction of the neutrino 
opacity at high densities due to rapid nucleon spin fluctuations p^- p!6| , p6| was explored. A 
suppression factor of 50% was assumed, which is the maximum reduction still allowed by 
consistency between the SN 1987A neutrino detections and the calculated neutrino signals 
for axionless cooling models Finally, the sensitivity of the derived axion mass limits 



on the mass of the cooling neutron star was tested by performing cooling calculations not 
only for our standard protoneutron star model with Mns,b = 1.53 M©, but also for models 
with baryonic masses Mng^b = 1.30 M©, 1.40 M©, 1.65 M©, and 1.75 Mq. Like most of the 
other comparative computations, the simulations with different neutron star masses were 
done with the reference values gap = 1.5 x 10~^° and x = 0. 



B. Protoneutron Star Evolution without and with Axion Emission 



For sufficiently large axion-nucleon coupling the cooling of protoneutron stars is signifi- 
cantly affected by the production of freely escaping axions. We compare here the results for 
cooling calculations of the 1.53 Mq neutron star, once without axions {gap = gan = 0) and an- 
other time with the representative choice for the axion coupling constants of gap = 1.5 x 10~^° 
and X = gan/ gap = 0. The energy loss rate due to axion emission was prescribed according 
to Eq. (p!2|), i.e., suppression or saturation of the axion emission at high densities and tem- 
peratures was not taken into account. The temperature evolution for the axionless case is 
shown in Fig. ^, the corresponding information for the case with axions is given in Fig. |^. 

The initial temperature profile {t = 0) corresponds to the situation a few ten milliseconds 
after core bounce. The temperature shows a fiat hump between the (baryonic) mass coordi- 
nates of M(r) = N},{r)m ~ 0.3 Mq and ~ 0.9 Mq {Nf,{r) is the total baryon number inside 
radius r). The temperature in these intermediate layers of the protoneutron star is higher 
than near the center because of the heating caused by the forming supernova shock. More- 
over, as the deleptonization of the protoneutron star progresses, electron degeneracy energy 
which is released in the process e~ + p —>■ n + u^, is not completely radiated away from the 
surface of the star by the emission of all kinds of neutrinos. Instead, due to downscattering 
and multiple absorption and reemission of diffusing neutrinos, a part of this energy stays in 
the star and leads to heating of the gas. Therefore the temperature in the interior of the 
star rises during the first ~ 7 s of the evolution (Fig. |^) and the temperature peak advances 
inward to the center of the star, following the motion of the layer where most of the lepton 
(electron) loss occurs. At the time when the temperature maximum has reached the center 
of the star (after ~ 7 s), the chemical potential fi^^ = He + jJ^p — jJ-n of electron neutrinos has 
dropped to its final equilibrium value = that corresponds to deleptonized, neutronized 
conditions. Now the star begins to cool down essentially coherently by the continuing energy 
loss due to the emission of neutrino-antineutrino pairs created by thermal processes. After 
about 38 s the temperatures are below ~ 3 MeV in the whole star. 
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If axions are produced by nucleon-nucleon bremsstrahlung in significant amounts, two 
consequences follow for the cooling. On the one hand, the maximum temperature in the star 
reaches only about 35MeV at the mass coordinate M(r) ~ 0.7 Mq (Fig. ^). This has to be 
compared with the peak temperature of ~ 48 MeV realized at the center of the star for the 
case without axion emission (Fig. 0). On the other hand, the temperature starts to drop right 
from the beginning and already after 18 s the whole star has cooled down to a temperature 
of less than 3 MeV everywhere. Obviously, axions are very efficient in transporting away the 
heat produced after the conversion of degenerate electrons into neutrinos. 

In Fig. ^ the local energy emission rate due to nucleon-nucleon axion bremsstrahlung 
is plotted against the enclosed baryonic rest mass M(r) for different times. The axion 
emission rate per baryon peaks where pT^'^ has a local maximum and is largest during the 
first few seconds of the evolution. At later times the temperatures in the star have dropped 
appreciably and the energy loss in axions is greatly reduced. This is underlined by Fig. |^ 
which shows the axion luminosity as a function of time in comparison with the combined 
luminosities of Ve and z/g (denoted by "Lj^^"), the combined luminosities of z/^, z/^, Vr and z/^ 
(labeled by "Lj^^"), and the total luminosity for all neutrinos {Ly). The axion luminosity 
decreases by an order of magnitude within 5 s, a time after which Ly has dropped by only 
a factor of 3-4. Note that the luminosities in Fig. |^ include reductions due to gravitational 
redshift and time dilation for an observer at rest at infinity. Since neutrinos diffuse out 
through the star and decouple from the stellar matter near the surface of the protoneutron 
star, their luminosities are corrected for the gravitational redshift at the stellar surface which 
is typically 25% (of the energy measured at infinity). Axions, instead, leave the star from the 
deep interior where most of the axion production takes place (see Fig. ^ and their typical 
redshift is 40-50%. Time dilation (i.e., redshift of the inverse time interval or frequency) is 
accounted for by the same factors. Therefore the amount of energy transported away from 
the star by axions relative to the energy emitted in neutrinos is larger than suggested by 
the integration of the luminosities depicted in Fig. ||. 



C. Parameter Studies 

Increasing the axion-proton and axion- neutron couplings gap and Qan^ respectively, leads 
to a decrease of the energy radiated in neutrinos. Fig. ^ displays the total energies lost by 
the 1.53 M0 protoneutron star in axions {E^) and in neutrinos {E^) as functions of the cou- 
pling parameter gap- The nucleon-pair axion bremsstrahlung is again described by Eq. (p!2|). 
The symbols mark computed models and are connected by cubic spline interpolation. The 
different lines correspond to the values x = gan/dap = 0, ±0.5, ±1.0. The axion energy is 
larger for larger absolute values of x and a slight difference between the results for positive 
and negative x (with the values of E^ being a bit larger for positive x) reflects the asym- 
metry of the last term in Eq. (|12D against changes of the sign of gan- In Fig. ^ both E^ 
and E^ are the energies as measured by a locally inertial observer at rest at the surface of 
the protoneutron star. The axion energy is therefore corrected for the gravitational redshift 
between the layers of the axion production and the neutron star surface. 

Fig. 1^ shows the number of z/g absorption events predicted for the KII detector, A^kii, 
and Fig. p| the time interval ^kii within which 90% of these events are registered. Figs. 



and 10 give the corresponding information for the 1MB detector. Because the fraction of 



the gravitational binding energy emitted in neutrinos decreases when more energy is carried 
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away by axions (Fig. P), all quantities displayed in Figs. [^-[T0| exhibit a rapid reduction 
with increasing gap- This reduction is stronger for larger For an axion-proton coupling 
constant between 1 x 10~^° and 2 x 10~^°, the expected numbers of detector events and 
the detection times are reduced to about half of the values for the axionless case. While 
the calculated 11 neutrino events within ~ 15 s for the 1.53 protoneutron star are in 
reasonable agreement with the KII measurement of SN 1987A neutrinos (11 neutrinos in 
12.5 s), the ~ 6 predicted 1MB events within ~ 9 s are clearly on the low side of the detection 
rate of SN 1987A neutrinos in 1MB (8 neutrinos in 5.6 s). This trend of our results is also 
present for supernova models of other groups and reflects the marginal consistency observed 
between the 1MB and KII neutrino data: For the z/g emission characteristics (spectrum and 
luminosity) deduced from the KII measurement, one would expect a much smaller number 
of 1MB events, or, inversely, Ve emission that has caused the 8 events in 1MB should have 
produced a much larger signal in the KII detector (see, e.g., ||36[| ) . 



If axion emission is absent, the expected event numbers A^kii and A^imb as well as the 
detection times txii and ^imb are slowly increasing with the mass of the protoneutron star 



25|. Between Mns,b = 1-30 Mq and Mns^b = 1-75 Mq, iYxn and Ajmb roughly double, 



while txii increases by a factor of 2.5 and ^imb by about 60%, as can be seen in Table 1 



of |25|. When axions are produced in the stars by nucleon-pair bremsstrahlung [Eq. (|T2|)] 
with axion-proton coupling gap = 1.5 x 10~^° and axion-neutron coupling gan = (i.e., 
X = 0), the results for different neutron star masses are obtained as shown in Fig. All 
four quantities reveal a very weak dependence on the baryonic mass of the neutron star. 
For the used values of the axion-nucleon coupling constants, we find A'kii — 5, Ajmb — 3, 
txii — 6 s, and tjMB — 3.8 s in all cases. Independent of the protoneutron star mass, these 
results are about half of the values found for the case of the 1.53 star without axion 
emission. The same inert behavior is also observed for the average energies of the electron 
antineutrinos captured in the two detectors: (ckii) — 23.7 MeV and (eiMB ) — 34.5 MeV. In 
more massive protoneutron stars the temperature and density become higher; this causes 
enhanced emission of axions so that a smaller fraction of the released gravitational binding 
energy ends up in neutrinos. The weak variation of the expected neutrino signals suggests 
that the additional gravitational energy release associated with a larger neutron star mass 
is essentially completely carried away by axions. 

Our study with varied parameters gap and x, based on the 1.53 protoneutron star 
model, was repeated for the case that the axion emission rate saturates at high densities 
and temperatures (see Sect. [Ill C] ). Employing Eq. ( P^ ) with Q*^^ from Eq. (|I2|) (with 
^ = 0.5) and F^^) from Eq. (|T|), the axion emission rate reaches its saturation level when 
PmTiV^ ^ 2W^(10MeV)-i. For an average interaction energy per nucleon of = 10 MeV 
we found that the corresponding reduction of the axion production increases the predicted 
number of neutrino events A^kii and Ajmb by roughly 40% compared to the number of 
expected events for the "naive" (unsaturated) case for which formally W oo. This trend 
turned out to be even somewhat stronger for the detection times txii and ^imb where the 
increase was about 60%. For W = 20 MeV the increase was ~ 20% in the event numbers 
and ~ 30% in the detection times, and for W = 5 MeV the event numbers and detection 
times rose by 60% and 80%, respectively. 

Taking into account the possible existence of a large number of negative pions at high 
densities in the simple way described in Sect. [11 B| and using the energy loss rate due to 
pion-axion conversion as given in Eq. (|H|), a reduction of the predicted detector response 
(A'kii, A'imb, ^kii and ^imb) to 75% of the values for the axionless case occurred for (jaN — 
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0.15 X 10"^°. A reduction by 50% was seen when QaN = (0.30 ... 0.40) x 10"^°, and only 
25% of the neutrino events were measured in roughly four times shorter time when cjaN was 
between 0.75 x 10~^° and 1.0 x 10~^°. Of course, the neutrino emission characteristics of the 
simulated protoneutron stars vary with time. Therefore, the detection rates in KII and 1MB 
are not constant, and, correspondingly, the values of the axion-pion coupling (jaN which lead 
to a certain reduction are somewhat different for the two experiments and also for the event 
numbers and detection times. 

Our results exhibit an extremely weak dependence on the choice of the nucleon degener- 
acy parameter ^ of Eq. (|n|). Changing ^ from (degenerate nucleons) to 1.0845 (nondegen- 
erate nucleons) leads to an increase of A^kii by ~ 0.5 events and to a detection time txii that 
is ~ 0.5 s longer. For the 1MB detector the corresponding numbers are ~ 0.25 events and 
~ 0.25 s, respectively. A larger value for ^ reduces the energy-loss rate by axion emission 
[Eq. (pr2[)], but the particular choice of is irrelevant at the level of accuracy implied by the 
sparseness of the neutrino data of SN 1987A. 

Replacing the hyperonic EOS B by the n-p EOS A in our 1.53 Mq protoneutron star 
model leads to essentially no change in the predicted neutrino signals. Hyperonization plays 
an important role only in stars with baryonic masses ^ 1.70 Mq where the density in a 
larger part of the star is high enough to favor the production of hadronic states other than 



n and p (see p2[ and compare also models S4BH_0 and S4AH_0 of Table 1 in [p5|] ). 



A suppression of the neutrino opacities relative to their "standard" values, e.g., by the 
many-body effects discussed in Sect. ^11 C[ can have an important impact on the predicted 
neutrino signal p5| and thus on the allowed range of axion-nucleon couplings deduced from 
the SN 1987A neutrino data. As a test, we reduced the axial vector contributions of neutral- 
and charged-current processes by about 50% (i.e., we chose a = 0.5, see [^]) in the case 



of the 1.53 Mq protoneutron star with axion emission prescribed according to Eq. (p!2| ) and 
Qap = 1.5 X 10"^'' and x = 0. The characteristics of the detector signals change in the same 
way as found for the axionless case in |2^ . A 50% reduction of the neutrino opacities leads 
to an increase of the predicted numbers of neutrino events iVxn and A^imb by ~ 2.2, but the 
detection times txii and tiMB drop by about 25%. The mean energies of Ue registered by 
the detectors rise by roughly 2 MeV. Since these changes are not dramatic, it is clear that 
even a 50% reduction of the neutrino opacity in the protoneutron star will not completely 
alter conclusions on the axion production drawn from a comparison of theoretical neutrino 
signals with the SN 1987A neutrino data. Although the increase of the number of detector 
events in case of a lower neutrino opacity can somehow compensate for the effects of axion 
emission, the detection times shrink both by a reduction of the neutrino cross sections and 
by the additional axion cooling of the star. Therefore we conclude that a value of the 
neutrino opacity that is lower than the "standard" one can hardly lead to a restoration of 
the compatibility between the SN 1987A neutrino signal and the signal predicted in case of 
strong axion emission. 



D. Excluded axion couplings and axion mass limits 

The cooling sequences of the 1.53 Mq protoneutron star model for varied axion-nucleon 
coupling parameters gap and x = g an /gap are used to construct exclusion curves in the x- 
(/ap-space. The sensitivity against changes of the axion couplings was found to be somewhat 
different for predicted event numbers and detection time scales. Nevertheless, it is safe to 
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claim that consistency between the theoretical neutrino signal and the SN 1987A neutrino 
data requires that A^kii,imb ^ i^Kii,iMB and tKii,iMB ^ |^fcii,iMB when A^kii.imb and tKujMB 
are the expected signal parameters for the axionless "standard" 1.53 Mq protoneutron star 
model. According to this criterion, the bell-shaped curves in Fig. ^ separate allowed from 
forbidden regions. The lower curve corresponds to the case where the energy loss in axions is 
described by Eq. (|12|), the upper curve represents the case where saturation of axion emission 
is included [Eq. (p4D] with an average interaction energy per nucleon of = 10 MeV. The 
values of the upper curve can be obtained from those of the lower curve by scaling with a 
factor of about 1.9. 

On the hatched sides above the curves the axion-nucleon couplings are so large that more 
than about half of the gravitational binding energy of the neutron star is emitted in axions 
and only less than half in neutrinos (compare Fig. There is only a slight asymmetry 
between positive and negative values of x which is caused by the minor asymmetry of the 
third term of Q^^ in Eq. (^) against changes of the sign of Obviously, the contribution 
of the asymmetric term oc QanQap is rather small. For this reason, the lower curve of Fig. ^ 
can be pretty well fitted by the relation g^^ + 2g1^ ~ (1.5 x 10~^°)^ which means that 
the excluded values of the axion-nucleon coupling constants lie outside of an ellipse with 
semiaxes 1.5 x 10"^° and (1.5/^2) x 10"^° in the fi'ap-S'an-plane. 

The deduced limits for allowed values of the coupling constants gap and x = gan/dap 
might be sensitive to the protoneutron star mass. This was tested by varying the baryonic 
mass of the protoneutron star for g^p = 1.5 X 10"^° and x = 0. In case of the 1.53 Mq model 
this pair of values {x,gap) lies on the lower exclusion curve of Fig. |T2|. The results of such a 
study are displayed in Fig. [Tl|. The observed inertia of the predicted neutrino signal against 
changes of the mass of the axion-emitting star suggests that the exclusion curves should not 
depend strongly on the (unknown) exact mass of the neutron star born in SN 1987A. 

Taking into account the possible existence of a large number of negative thermal pions 
in the protoneutron star in the simple way described in Sect. [11 B[ we find that the rapid 
energy loss by axions emitted in pion-axion conversion processes [Eq. (p^j excludes couplings 
of roughly cjaN ~ (0.30 ... 0.40) x 10~^°. This assumes a perturbative treatment of the 
emission processes and neglects possible saturation effects. Since for typical conditions 
in the supernova core the perturbative energy emission rate by pion-axion conversion is 
approximately 10-20 times larger than the lowest-order energy-loss rate for the nucleon- 
nucleon axion bremsstrahlung, the bound on the coupling constant cjaN would be 3-5 times 
more stringent than the limit on gap in the perturbative case. 

With these limits on the axion coupling parameters gap, gan (or, equivalently, x), and 
(jaN, upper limits on the axion mass can be derived from Eqs. (^ and (|I3]) when use is 
made of the experimental values for the dimensionless couplings Cp and C„ (see Sect. [11 A] ) 



collected in Fig. |I], and Eq. (|I|) is employed to relate the PQ scale (or axion decay constant) 
fa with the axion mass rUa- From Eqs. (^ and (|I|) one finds 

ma < 0^ ^25) 



10-3 eV ~ Cp 10-10 



and from Eqs. (|T5|) and ([^) one obtains 



nia ^ 0.66 (1 ^ x'^ x\ g'^§'^ 



10-3eV a \2 2 3 / 10-10 



(26) 
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Here g^p^ and g^§^ are the maximum values of the axion couphngs allowed by the SN 1987A 
neutrino detections, and Cp = fi'ap/["^p/(/pQ/^)] and x correspond to the pair of coordinate 
values of a particular chosen point in Fig. |l|. Note that g^p^ as a function of x describes the 



(lower) exclusion curve in Fig. 112 



The upper bounds of the axion mass deduced for the points of Fig. |l] by this procedure 
are plotted in Fig. |13|. Above a line connecting the point {gan/gap,^a/[^0~^ eV]) = (—1, 1) 
with the point (3, 5), Fig. [1^ shows the results for the "standard" case without pions, whilst 
below this line the dots represent the mass limits for the case with pion-axion conversion. In 
both cases saturation effects have been neglected. The latter case yields axion mass limits 
that are stronger by a factor of almost four. As was mentioned in Sect. [Ill C| , the inclusion 
of saturation effects is expected to lead to limits that are less sensitive to the presence of 
pions. 



V. SUMMARY 

We have re-examined the stringent limit on the axion mass inferred from neutrino emis- 
sion by SN 1987A, in the light of a possible suppression of axion emission by the many-body 
effects of nucleon spin fluctuations and additional emission processes involving pions, taking 
into account the latest determinations of axial-vector current couplings to nucleons. The 
suppression of axion emission due to many-body effects degrades previous limits by a factor 
of about 2. Emission processes involving thermal pions can strengthen the limits by a factor 
of 3-4, if saturation effects on the nucleon spin fluctuations are neglected, whereas inclusion 
of such effects tends to make the limits less sensitive to pion abundances. The resulting axion 
mass limit depends upon its precise couplings, ranging from 0.5 x 10~^eV to 6 x 10"'^ eV. 
Our results are consistent with previous limits []TU[, though more precisely stated. Fig. W2 



shows our limit on the axion-proton coupling as a function of the ratio of the axion-neutron 
to axion-proton couplings, and Fig. ^ shows the mass limit as a function of this same ratio. 

To test the stability of the deduced axion mass limits against uncertainties of the stel- 
lar models and of the description of the input physics, we performed a large number of 
comparative computations. It turned out that neither a change of the nucleon degeneracy 



parameter ^ of Eq. (pAf ) nor a reduction of the neutrino opacity to 50% of its standard value 



lead to major changes of our conclusions. There are differences of our work compared to 



previous work by Burrows et al. ||T0[ concerning the protoneutron star modeling, namely we 
used different equations of state and did not take into account accretion of matter onto the 
nascent neutron star. 

Two different equations of state were used in the presented models, an n-p equation 
of state and an equation of state where hyperons occur as additional hadronic degrees of 
freedom at densities beyond about two times nuclear matter density . For neutron stars 
with baryonic masses ^ 1.70 M0 hyperons are abundant only in a relatively small central 
core region of the forming neutron star and therefore do not dramatically change the cooling 
and neutrino emission of the star. Similarly, the axion emission was found to be affected 
only at a minor level by the presence of hyperons. 

The amount of material possibly accreted onto the protoneutron star during the first 
moments of the supernova explosion and the corresponding accretion rate as a function of 
time are rather uncertain and must depend on the structure of the progenitor star and on 
the details of the still incompletely understood explosion mechanism of Type II supernovae. 
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Instead of introducing new free parameters to model accretion, we tried to account for the 
unknown amount of accreted matter by performing cooling simulations for protoneutron 
star models with baryonic masses between 1.3 Mq and 1.75 Mq. Interestingly, the higher 
temperatures in more massive stars lead to enhanced axion emission that carries away a 
larger fraction of the gravitational binding energy that is released during the cooling of the 
star. Therefore the predicted neutrino signals in the KII and 1MB detectors and thus the 
axion mass limits deduced from the SN 1987A neutrino data should be rather insensitive 
to the exact mass of the protoneutron star formed in SN 1987A. Including accretion as 
in Burrows et al. |]IU| raises the total neutrino luminosity above the contribution from the 



core and leads to additional events in the detectors which are essentially unaffected by the 
emission of axions. Therefore the event numbers Nym and Niub in reference were found 
to drop somewhat less strongly with increasing axion-nucleon coupling than they do in case 
of our models. 
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FIG. 1. Scatter plot|^c^l^^ ixicJi^/otHr?" Coupling an^the ratio of the axion-ncutron to ax- 
ion-proton coupling for different axion models allowing for the uncertainties in Aq; from the bottom 
left to the top right, /? = 0, /? = 27°, KSVZ, P = 54°, and /? = 81°. Note that the axion-neutron 
coupling is much smaller than the axion-proton coupling for the KVSZ axion and for the DFSZ 
axion when /? ~ 45°. 



20 




FIG. 2. Temperature evolution of a 1.53 protoneutron star without axion emission. The 
temperature profiles are plotted against the enclosed baryonic rest mass M{r) = Njj{r)m (in solar 
masses; m is the common nucleon mass) for different times from the start of the computation 
shortly after the formation of the protoneutron star. For time t = the curve is labeled, the 
other times are listed according to the order of the corresponding curves (from top to bottom at 
M = IMq). The interior of the star heats up first due to the conversion of lepton degeneracy 
energy into thermal energy, but finally cools. 
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FIG. 3. Temperature profiles of a 1.53 protoneutron star vs. enclosed baryonic mass for 
different times. Axion production via nucleon-pair bremsstrahlung is included with axion coupling 
constants gap = 1.5 x 10~^'^ and x = gan/dap = 0. The star cools much faster than in Fig. ^. 
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FIG. 4. Local energy production rate per baryon due to nucleon-nucleon bremsstrahlung 
emission of axions [Eq. ([l^)] with coupling constants gap = 1.5 x 10^^*^ and x = gan/dap = 0. 
The profiles are plotted against the baryonic mass coordinate of a 1.53 Mq protoneutron star for 
different times after the start of the simulation (the times are listed according to the order of the 
curves at M = 0.5 Mq). The energy production per baryon roughly peaks where the product pT^-^ 
reaches a maximum. 

23 




FIG. 5. Electron neutrino plus electron antineutrino luminosity (labeled by L,^^), sum of 
all luminosities of heavy lepton neutrinos (labeled by Li,^), total neutrino luminosity Ljy, and 
axion luminosity as functions of time during the cooling of a 1.53 protoneutron star for 
axion-proton coupling gap = 1.5 x 10"^'^ and axion-neutron coupling gan = 0. The luminosities are 
redshifted as measured by an observer at rest at infinity. The axion luminosity drops by an order 
of magnitude within 5 s, whereas the neutrino emission decreases less rapidly. 
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FIG. 6. Energy loss of a 1.53 Mq protoneutron star in neutrinos {E^) and in axions {E'^) for 
different combinations of values of the axion-proton coupling gap and axion-nGutron coupling Qan 
m Eq. (H). The energies are given as measured by a locally inertial observer at the surface of 
the neutron star. The symbols mark computed models, the interpolation is done by cubic splines. 
The curves through the asterisks correspond to the case x = gan/dap = 0, the circles to the cases 
X = ±0.5, and the squares to x = ±1.0. Filled symbols mark positive, open symbols negative 
values of x. 
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FIG. 7. Prediced number of Ve absorption events in the KII detector for a 1.53 Mq protoneu- 
tron star as a function of the axion-proton coupling gap. The symbols correspond to computed 
models and the different curves represent different values of x = gan/Oap (see Fig. ^). 

26 




FIG. 8. Time intervals txii where 90% of the predicted capture events in the KII detector 
happen for a protoneutron star with a baryonic mass of 1.53 and for different combinations of 
values of gap and x = gan/Oap- The meaning of the symbols is the same as in Fig. ^. 
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FIG. 9. Same as Fig. |, but for the 1MB detector. 
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FIG. 11. Numbers of predicted absorption events in the KII and 1MB detector, Nym 
and A^'iMB) respectively, and corresponding time intervals tKii and tiMB where 90% of these 
events happen. Five protoneutron star models with different baryonic masses, M^s^b = 1.30 Mq, 
1.40 Mq, 1.53 Mq, 1.65 Mq, and 1.75 Mq were evolved with axion emission by nucleon-nucleon 
bremsstrahlung [Eq. (|l^)] for gap = 1.5 x lO^^'' and x = gan/Oap = 0. The expected measurements 
for event numbers and detection times vary extremely weakly with Mns,h- 
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FIG. 12. SN 1987A limit to the axion proton coupling as a function of the ratio of the 
axion-neutron to axion-proton couplings, with and without saturation due to many-body effects. 
Saturation (for an average interaction energy per nucleon of = 10 MeV) degrades the limit by 
a factor of about 1.9 independent of the ratio of the axion-neutron to axion proton couplings. The 
combinations of coupling parameter values above the bell-shaped curves are excluded. 
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FIG. 13. Sn(iA^ IguiftJ axK)?^a^lss^Sl^^S^on of the ratio of axion- neutron to 
axion-proton couplings for different axion models, with (lower set) and without (upper set) pion 
processes; from bottom left to top right, /3 = 0, /? = 27°, KSVZ, (3 = 54° and (3 = 81°. Note that 
saturation due to many-body effects has not been included (it degrades the limit by a factor of 
about 1.9). 
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